clear;clc
% f=@(x)x^4-4*x^3-6*x^2-16*x+4;
% x0=fzero(f,[-1,4])
% x1=fminbnd(f,-1,4)
% f(x1)
f=@(x)x(1)^2+25*x(2)^2-2*x(1)+50*x(2)+25;
x0=fminsearch(f,[2 3])
f(x0)